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A.  Theory  for  Strong  Electric  Correlations.  We  have 
developed  a  theory  for  the  probability  distribution  function  of 
the  electric  field  in  a  highly  turbulent  environment.  This  makes 
contact  with  our  previous  experiments,  in  which  we  found  a 
characteristic  distribution  function,  exp(-E  )  in  a  range  of 
different  beam-plasma  strengths.  The  theory,  as  aetailed  in 
Appendix  A,  extends  previous  work  which  took  into  account 
collective  effects  in  the  plasma  surrounding  the  ion.  uur  work 
takes  into  account  the  possibility  of  strong  wave  correlations 
overlengths  exceeding  the  Debye  length,  k^.  We  find  that  to 
explain  the  lOOkV/cm  fields  we  have  observed  we  must  invoke 


correlations  which  are  quite  strong  over  dimensions  of  at  least 


10\D.  This  is  a  new  result  with  possible  powerful 


consequences.  It  makes  connection  with  the  general  body  of 
soliton  theory,  which  envisions  the  final  stage  of  plasma 
turbulence  as  compressed  packets  of  electric  field  on  the  scale 
of  10-100  Debye  lengths,  and  containing  field  strengths 
comparable  to  the  local  thermal  pressure.  The  fact  that  we  have 
found  this  distribution  in  a  very  general  way  implies  that  it  may 
be  a  common  feature  of  strong  turbulence.  Further,  it  confirms 
our  previous  assertion  that  t no  specific  form  expt-E^)  implies  a 
one  dimensional  nature  of  tr.o  electric  field  spectrum.  If  the 
electric  fields  were  spheric  lily  symmetric,  for  example,  the 
distribution  function  would  E“exp(-E)  .  Tins  implies  a 
picture  of  one-dimensional  .  .--pi  isma  turau fence  whicn  persists 

for  the  entire  duration  of  - ...  r  «.•  a  r.  pul.-:  ,  at  Least  a 
microsecond. 


B.  Optical  Study  of  electric  fields  less  than  lOkV/cm.  We 
have  carried  out  some  of  the  studies  we  promised  in  our  last 
proposal,  of  the  distribution  function  for  relatively  weak 
electric  fields.  This  is  treated  in  detail  in  Appendix  B  and 
shows  that  we  have  successfully  developed  the  diagnostics  for 
this  sensitive  regime.  We  find  that  the  distribution  function 
remains  the  same  at  lower  electric  field  strengths.  We  are  now 
undertaking  to  measure  the  electric  field  distribution  as  a 
function  of  distance  along  the  propagation  axis,  radius  out  from 
that  axis,  and  time  during  the  pulse(z,  r  and  t) .  This  will  form 
the  principal  thesis  topic  of  a  graduate  student.  Ami  Dovrat. 

C.  Tne  anisotropic  nature  of  turbulent  microwave  emission. 

We  nave  carried  out  preliminary  studies  of  the  anisoptropy 

of  microwave  emission.  This  snows  that  there  is  unexpected  and 
quite  possibly  meaningful  dependence  of  the  polarization  of 
radiation  on  the  angle  of  emission.  We  hope  to  use  this  to  probe 
the  nature  of  the  emitting  entities  deep  within  the  plasma. 

Since  we  suspect  these  regions  are  quite  small,  on  the  range  of 
100  Debye  lengths  or  less,  radiation  remains  the  easiest 
diagnostic  of  such  features.  We  suspect  that  collapsea  or 
"caviton"  entities  produce  the  radiation,  and  are  generally 
oriented  along  the  magnetic  axis,  whicn  is  the  same  as  the  beam 
axis.  However,  the  polarization  pattern  does  not  correspond  to 
such  a  simple  one  dimensional  picture.  we  shall  try  to  interpret 
tile  data  in  terms  of  a  statistical  distribution  of  dipoles  witnin 
the  plasma.  This  is  being  undertaken  now  botn  ex  per imentnl ly  and 
theoretically. 
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Appendix  A 


Electric  Microfield  Distribution  In  Plasma 


With  Long-Range  Correlation 


Gregory  Benford  and  Xiaoling  Zhai 


October  2,  19S7 


Collective  and  individual  particle  correlations  affect  the  probability  di'ti 
bution  W(E)  of  the  electric  microfield  in  a  stationary,  turbulent  plasma  \\. 
extented  previous  work  to  include  long-range  correlations  over  lengths  »  \ j , 
the  Debye  length.  The  characteristic  distribution  U'(  E)  a  exp(  -  1 

emerges,  with  d  the  dimensionality  of  the  electric  field.  I  he  mean  ■  I  J  ■  is 


proportional  to  the  square  of  the  correlated  particle  density  In  descr 


lie'  a  re.  .  n I 


strongly  turbulent  experiment  with  <  E~  t*5H‘/crn  requires  .  ,  ,rre|ut  mn 

over  scales  >  lOA/y.  Comparison  with  observed  microwave  emission  implies  that 
more  than  one  percent  of  the  plasma  volume  experiences  strung  field  regions 


V  ./vt  O  V.v.  c.  f.V.  C.  V.  *  c. 


1.  Introduction 


pi iwwvtm’mmv.m 


r  s  w  v  v.rfvc  w  v.  v 


V.  V.  V-  V.  V.  /.  V,/. 

V. 


Many  experiments  have  measured  the  electric  field  distribution  inside  tur¬ 
bulent  plasma.  (Klein  and  Kunze  1973  ;  Gallagher  and  Lavine  1971  ;  Hamada 
1970  ;  Antonov  et  al  1970  ;  Matt  and  Scott  1972).  A  recent  measurement  (Lev- 
ron,  Benford  and  Tzach,  1987)  found  a  probability  distribution  of  the  electric 
Langmuir  field  in  a  strongly  driven  beam-plasma  environment, 


VV(E)  a  exp(—  - 


■) 


<  £2  > 

where  <  f2  >  is  about  (85kV/cm)2. 

Theoretical  investigation  (Ecker  and  Spatschek  1973;  Hooper  1967;  Ecker 
1971)  of  a  stationary  equilibrium,  homogeneous  plasma  model  gave  some  phys¬ 
ical  explanation  of  the  electric  field  distribution  in  terms  of  screening.  All  yield 
an  exponential  distribution.  G.  H.  Ecker  et  al  (1973)  give 

lT(E)ocexp(  — 11) 

Eh 

where  £’//  is  the  Holtsmark  field  (Holtsmark  1919).  It  is  about  0.16k V/cni  in 
the  experiment  of  D.  Levron  et  al,  about  10-3  times  smaller  than  observed. 
Plainly  collective  effects  and  possibly  long-range  correlation  must  conn*  into 
play  to  explain  this  striking  difference.  Generaly  for  turbulent  elect ri*'  field  of 
dimensionality  d,  \V{E)  oc  Ed-1exp  (-E2).  The  I.evronet  al  experiment  implies 
d=l.  with  such  strong  antisotropy  apparently  arising  from  the  beam-plasma 
instablity. 

Ifer<’  we  develop  a  theory  for  U  (E)  which  decomposes  tin-  e|,  ctric  mn  ioli.  Id 


distribution  into  three  parts:  individual  particle,  collective  and  long-range  corre¬ 
lations,  We  give  each  an  equal  footing,  following  the  strong  wave  field  approach 
of  Ecker  and  Spatschek  with  correlations  added. 


2.  Basic  Equations 


Let  us  start  from  the  Maxwell  equations,  with  no  net  current  j  =  0.  B  =  0 


and  a  homogeneous  time-independent  plasma. 


V  •  E  =  £(r  —  r, )  V  x  E  =  0 


V  x  B  =  0 


V  B  =  0 


We  can  define  a  vector  potential  A  with  Fourier  components 


A  =  c\j Y  A**‘ko  r 

k 


With  ko  =  j-,  k  =  |k|  ,  and  since  VxE  =  0.  E=V  v.  with 


--'VrETf,kr 

k 


Using  equation  (4)  in  equation  (l),we  can  get: 


ie  4*  ^  lkr 


it:  /  h  <t 

T\l~ ^ 


The  Hamiltonian  of  our  system  in  the  random  phase  approximation  is 


h  =  V^,iv 


El1,'  -r-T  y,  y'  y^  J_ 

•’,n  *"  r  2^  r- 


,k  r,-r,  i 


We  are  free  to  choose  .4k,  so  we  choose  .4k  =  0  at  the  second  step  of  equation 
(6).  The  Hamiltonian  has  the  canonical  conjugate  variables,  r, ,  p, ;  .4k,  E jc. 
The  Hamiltonian  equations  are: 


4 


<)H 


(7) 


Consider  the  ensemble  distribution  function: 

i=N 

Fr=  I]  «(r’  -r.(t))6(p‘  -P.(0)Il6(Ak  -4)6(£’k  ~  4(0)  (8) 

«'=i  k 

Where  r‘,  pl,  .4^,  E k  are  the  coordinates  in  T  space,  and  r, .  Pi,  -4k.  Fk 
the  coordinates  in  /j  space.  We  define 


Hr  -  r,-(<))«(p  -  Pi(0)  n(^k  ~  4W£k  -  4)  (9) 

i  k 

Here  the  sum  is  over  all  particles.  Now  consider  the  time  derivative  of  F1A, 


,dFi 


=  (^) 

dt  'r.p,(>tk  £k)  '  <9r,  'r.p.p.  iaK  fk,  ,4^  ■ 


+  E< 
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an 

9Ek 


an 


Vp.r,  (.iK 

,£kH4k,£] 

dFxl> 

<9.4 1 

^r.p.r,. 

p.  M*  £*) 

dFlu 

dE{ 

^r.p.r,. 

p,.(.4k  .£*) 

'~.T , 

u 

bv  &L 
-  op, 

.  P.  by,  - 

:  10) 


an 


it  in  equation  (6)  then  yields 


/■It 


<Tif>  - 

k 


«k  r. 


\ 


jjS 

w 

ll 


+  Kko  •  loVMlfk  +  »>  ^(k+1)  r'(d-£r 


-  rDylEoo  p. 


,  a^i 


mV  V 


+  ^Epoio)4^k+1-](^) 

Following  Ecker  et  al  (1973)  we  use  the  collective  coordinate 


«k  =  Vw 


Instead  of  putting  all  electron  motions  into  the  collective  coordinates,  we 
choose  a  value  of  k  such  that  the  collective  modes  only  enter  with  wave  vector 
less  than  k,-.  Clearly  kc  >  r^1,  and  for  undamped  modes  we  require  k  >  kc,  with 
kc  lying  between  r^1  and  Xp1.  Previous  calculations  took  kc Xq  «  1,  with  A £> 
the  Debye  length.  Here  we  explicitly  expand  the  formula  to  include  correlations 
for  k  <  Xp1  (Fig  1).  Define 


=  “p("S  E 


i,k,fc>fcc 


Make  the  transformation: 


P.  =  P<  - 


for  k  >  kr. 


Pi  =  s  p,5 

■4  =  5-1 4s 

(13) 

E  k°4eikr‘ 

(H) 

h\  =  £Y.  +  i\ 


.V  .V  A  jy.V.  .  '  .A .  »  ,VA  _■»  A .  .\  A  A  .V  A  A  A  A 


(16) 


for  k  <  ke,  Ajj  =  t 

3.  The  Distribution  Function 

We  break  the  Hamiltonian  equation  (6),  into  contributions  from  the  two 
salient  regions  of  k  space: 

*  -  Eg+^EE  E 

+  1/2  E  <£k£'-k  +  yAkA-k)  <17> 

k.i<tc 

The  plasma  wave  dispersion  relation  is 


>  - f - e  '  =  U  tisj 

)"2  =  1  (I9) 

mV  '  m 

I 

In  what  follows  we  frequently  assume  a  cold  plasma,  ((I'Pj/mw)‘)  1  Now 

we  approximate  the  complex  terms  of  equation  (IT).  Since 


2|°~T  y  y  y  _L  .kir.-r,)  _  _Ly  y  f _ 

v  \fi  ~  2V  ^  ^  |r 

1  }  k  k<kc  1  id>1 

Using  this  to  equation  (11),  through  some  algebra  we  get. 


f-*,|r,-r,| 


r,  -  r;  | 


dF,^  p.  ,9Fitl  v'  i  ^ 

— +  ™  y '  E.  rA-m 

K.k  <kc 


k  Tf/TT  + 


v  J  P  J  A 


di-dp  J]  dAxdE\  —  E  '  ;  -  0 


We  can  illuminate  this  complicated  result  by  using  the  1  space  density  Pi 


the  first  order  distribution  function  f,,.  and  the  second  order  distribution  funr 


tion  fij. 


LruHi'L 


'M 


Pd  =  Pd( ri.r2 . r.v.Pi.Ps . p.v, .  .4kj . t\  .  £k . 0(22) 

r  J  =  v 

/■<•  =  /  dr,  dp}  ]^[  f/.4k<//-.k  (23) 

7  j  =  i  k 

/  .  "=-v 

/■  j  c  =  J  F i  n  Pj  ii  Pd  n  dr„dp„  dAydE-y^  (24) 

U—  1 

where  N  is  the  number  of  electrons  in  the  system. 

Mulitplving  equation  (21)  by  Pd  ,  then  integrating  over  all  r„.  p„.  .1  k ,  £k. 
and  using  equation  (23)  and  (24),  yields 

f  +  &-LUJ*** 

-  Z  T^kl^+E^klr1^  (25) 

^  l  Kdtv  Kd4 

k  ,ifc<lfcc  K  k  K 

We  define  R  as  the  probability  function  of  collective  electric  field,  with  /, 
the  probability  function  for  the  individual  particle  electric  fields.  ;.c.  /,,  - 
/i/f.  fijt- —  f,jRR.  Since  and  /,Jt  are  separable,  we  ran  separate  equal  ion 


(25)  into  two  equations: 


‘>f>  |  P  Of,  c 
ill  m  Or 


i  /"  •  ■  <3  e  *'4i'  t'l  it 

-r  drdp—' — ■ - r— TC-(V  >  =  0 

J  Or  |r  —  r|  <tp 

< I R  OR  OR 

7 JT  ~  ^TAkoK  +  ^  'kool 


Kquation  (27)  can  he  formally  solved  as 


!{  =  n  h[  Ak  A-k'y  +  y  k  !  -  k  1 

k  *<*.. 


•  s  *»  s’  %"  %  %’  V  ■»'  ‘j 


I*  *.»  kH 


'#*1 


a 
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We  can  get  E(r)  from  equation  (5), 


e|'|=tSSt'‘i‘1"'1  c») 

'  k 

We  can  separate  E(r)  into  single-particle  and  collective  parts,  with  the  first 
term  counting  individual  particles: 


E(r)  =  E  ini  +  E  con 

.  4e  7T  t — >  r — '  k(; 

£  £  T 

1  k,t>tc 


}S£  „.k  (r-r.)  ,  .  /i 


T7  £  4k° 


r  (30) 


k.i:  <  trc 


Split  the  electric  field  Econ  into  two  parts: 


Econ  —  Ec0n  4-  Ecorr  — 


£  4k0e' 


E coll  represents  collective  screening  while  Ecorr  describes  correlation  over 
long  range,  k  <  kc.  Now  we  define  k V ( E )  as  a  probability  distribution  function 
of  the  electric  field,  with  P,n^  and  Rr0n  as  distribution  functions  of  Em(j  and 


W(E)  =  J  £(E  -  E,nd  -  Ecoil)RcollR,ndd[ . (32) 


Let  IV(q)  be  the  Fourier  transformtion  of  IV '(E) 

IV(q)  =  J  IV(E)  i  Ef/E 

=  J  f  _l<l  1  e.k  +  E,„.i  ]  Rron  Rlndd[....},n,iil[  (33) 

Since  R,„dd[. and  R.  .. •/'  by  assumption  vary  independently.  we 


can  write: 


V(q  )=/  E— -  /  <’'‘E . : 

J l  nd  '  ■  ii 


1  l  ,,,  •!  (J  I  1 1 


^  jr»  r*  v^jt*  ."*' \r*  vv  ^  vv  .^nj 


We  can  take  r  =  0  without  changing  the  electric  field  distribution.  Now 
consider  E ^  in  two  parts,  E^  =  E^rr  -f  E^11 ,  or  explicitly. 

-  ‘l.fE-'kr  l'v'"'  .-'it. 


4  = 

1=1  1= 1 

Here  :V,  represents  the  number  of  correlated  particles  and  ,Vnc  represents 
the  number  of  uncorrelated  particles.  The  spatial  dependence  of  correlation  is 
y(r,/L)  which  we  shall  choose  to  yield  finite  <  E~  >  in  a  correlation  scale  length 
L.  For  correlations  we  require  k  •  r,  =  2mi,  where  m  =  0;±1;±2;±3; . ,  so 


pcorr  _  ie  V  „tr* t  ~  ‘Uie  r *  f 

4  >*— vW». 


nry(r/L)r~dr 


Here  R <*.  =  ~ 1/3,  R0  is  a  very  small  but  finite  length,  and  nc  =  is  the 


density  of  correlated  particles. 


We  can  write  Eck°Tr  as. 


4orr  =  ?Vyn‘GW 


/■tto 

k)  =  4* 

Jrb 


y(-)r2dr 


E(r  )=yn'ie  Y,  fG{k)+\V  E  Ek* 


k.Ktc 

It  is  easy  to  see  that  =  —  Ejj^  and 


k  k<ke 


Ek"+wr  =  Vr  E  =  E  2*'<4> 

k  Kt,  k  *<t,. 

Where  E\Llt<k  includes  Ej{t>,,  and  Ej(orr.  now  nse 

Rro„+r,,rr=  n  /k(,lk,rk^  +  /;kE-k) 
k  k<k. 


&* 


■ 


VVrcoi/+conr  ( q)  is  the  Fourier  component  of  lVcon+corr  (E),  so 


"(E)  =  / 


e-qCjWcoll+COrr(ci)dq 


3  ^2 

WW.rr(E)  =  5nKHjcikdk)WeXp(-^fCikdk] 

Take  <  E2  >  =  32 it  f  C'lkdk  ,  then 

Wcoii+corr(E)  =  n/27/8^3/2  —  E2  >3/2  e2/J<Ej> 

—  3E2 

We  can  see  the  distribution  of  electric  field  is  proportional  to  e  2<e2>  %  where 
E  includes  Ecoi;  and  E<;0rf. 


E(r)  =  Ecorr  +  Ecoll  =  ynce  £  JLG(k)+yJy  ]T  4k°  (38) 

k.fc  <jfcc  ki<tc 

From  equation  (38),  we  can  derive 


EeoifE corr  =  (y)3/7nce  £  G(k)/*Ek 

k,*<tc 


-  ie  Mjr^r 

Ek  =  j\jy  L, 


i  =  1\nc  ;=.V„ 


=  £  E  E  E  e'(k  r'_'  r,)(ko  •  lo) 

k,k<kc  1  =  1  >  =  1 

is  the  number  of  uncorrelated  particles,  a  large  number  e-,flr' 

vanshes,  because  ik  ■  r,  is  a  random  phase.  If  i=j  and  k  =  —I,  then  we  have 
e-lkr.+lr,)  _  j  so 


E;oll  = 


1 1 


ELr 


Later  we  will  see  that 


<  El 


"2  E  £  t  ■  7G(fc)G(/) 

k,*<icl,;<tc 

>  is  a  large  number,  so  we  can  neglect  Eco;j  • 


Ecorr  and  find  E^ol(+oorr  as  E2^,  +  E2orr  Now  we  have,  for  a  one-dimensional 


field, 


Wcoll+c  orr(E)  = 


E;.„+Eg.rr 

g  2/3  <£*> 


<  EJ  > 


8  7T3/2  <  £2  >3/2 

=  ^  Lk<k  ^k^-k<ik  =<  e^"  > + <  E^rr  >  +6’ 


Here  the  residue  is 


" = +tAt''A  m 

Because  of  the  random  phase,  both  E^11  E™£  and  E^rr  E™£  become  negligible 
compared  with  <  2?2o((  >  and  <  £2orr  >,  so  l*=0  The  distribution  function  of 
electric  field  is  generally. 


Wcoll+cor  r(E)  = 


27 


1 


(•2  4.  F  2 


<E,Ol  (>  +  <W> 


8  7r3/2(<  £*,,  >  +  <  £20Pr  » 

We  are  interested  in  the  case  where  (£2orr)  >•  (E;olt).  so  U',.„a+rorr t E)  as 
Wrorr.  If  E,^,.,.  is  much  stronger  then  Eco|j,  the  distribution  funtion  is  nearly 
the  same  as  the  correlation  electric  field  distribution.  This  is  the  distribution 
function  for  one  dimension.  For  d  dimensions,  we  should  have 


W 


rd-l. 


For  three  dimensions  the  distribution  function  becomes 


U'  x  £V 


a  form  often  seen  in  weakly  turbulent  environments  (Klein  and  Kunze.  l!)7.'i). 


Examining  <  E;oll+corr  >, 


<  EL,  > 


4  ir  /  1 

/  vncpr A 

ik<kc 


4tr  e2  4jt 


<  Crr  > 


L  ?«*>** 


(41) 

(42) 


The  spatial  dependence  of  correlations,  j/(r),  allows  finite  <  E;orr  >  only  if 
y(r)  =  ( L/r )3'2.  This  is  an  artifact  of  our  method;  in  general.  y(r)  must  be 
rapidly  declining.  We  find 


G(k)  =  4jt  fR°°  (L/r)3t~r2dr  =  6ttL3,‘(V1/2  -  V01^2' 
Jro 


<  £-2  >  _  nF:4iL3(4T|34,  ,t  ,  to  ol0 

^  <•  nrr  ^ 


.1/2 


2t2  9  M1  +  v  2  l'!/2  J 


Here  V'o  =  Rq;  in  the  limit  V  -  x 


.  2  .  128  2  r,f3i 

*■-  kr/jrr  > —  ^  uce  L  k? 


Since  <  E;  u  >=  lfi;r2nn,.e2.  we  find 


/■oil 


,2 ,3, 


^  ^,/ro/i4-'*orr  ^  ^  r  ^  L  A “f*  f  7lrirk,- 
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(43) 

(44) 


(45) 


(46) 


He  re  —  A^n^/1  is  tin*  density  of  particles  which  are  uncorrelated. 

4.  Comparison  With  Experiment 

To  compare  with  observations  we  use  the  the  data  of  the  I).  I.evron  cl  al 
( 1987).  Their  density  of  part  icl*-s  is  about  10l3cm_  i.  mean  electron  temperat  ure 


is  about  lOeV.  and  the  Debve  lenet  h  A/j  ss  7.4  x  10  Vm.  so  l  /i  s;S  IxlllVm  1 


The  interparticle  separation  r0  is  about  1.3  x  10  5cm.  so  ko  —  jr  «  5 7 X' /j>  The 
Holtsmark  field  is  Eh  =  =  0.16kV/cm.  We  can  write,  with  t  =  ^-ru, 

ro  * 

<  £2  >=  (47) 

o  n 

Levron  e<  a/  found  <  £2  >0jaKs  (85kV/cm)2,  about  3  x  105  times  greater 
than  E'h.  To  fit  experiment,  we  need  ncL3)k  1.3  x  104.  Let  us  see  what 
this  implies.  Define  Xd  as  the  number  of  particles  inside  the  Debye  sphere 
( ,\’d  a*  1.5  x  10s).  .\'r  =  ^-nrL3  is  the  number  of  particles  inside  the  correlated 


sphere.  Then 


0.25 


I'  —  * 

Nd  n  ’  Ad  [t(^)2]  1/3 


(48) 


If  we  use  kr  —  ko,  k  is  about  0.02.  Even  if  we  choose  7^  =s  1 .  for  complete 
conhence  in  the  region.  L  is  about  Ad-  mor  possibly  taking  ^  =  0.1.  L  is  about 
5 Ad  Thus  long-range  correlations  are  necessary. 

We  can  use  further  constraints  set  by  the  observed  microwave  power  from 
the  same  experiment.  Using  the  electric  field  distribution,  we  can  calculate  the 
power  of  the  microwave  radiation  of  the  plasma 


P,,tai  =  X,  J  J  \V(E)dE^dQ 

Where  \V(E)  —  *J  .  is  the  electric  field  distribution  (one  dimension),  with 
u~  —<  E~  >  X,  is  the  number  of  cayitons.  We  take  the  power  emitted  by  a 
single  soliton  correlated  over  a  scale  L.  There  X,  such  emitters,  which  we 


represent  as  simple  dipoles  oscillators  of  dipole  length  L.  I  hen 


Ptotai  =  o~6n- X,ffTcEH~[nrL  )  {—  )k 

n 

I  he  experimentally  observed  result  (Barauga.  Benford.  Tzach.  and  Kato.  1085) 
is  about  P„b,  —  1 05 f\  Watt  Since  tjc  a:  5  x  II)-*3  Watt,  we  found 


(—)\nL'fS.k  a;  3  x  1U25P5 
n 


lake  nZ.3  «  5  x  lO4!^-)3.  plus  the  observed  u*  =  C)4/‘.iE]j^-(nrL3)k  = 

(85 H'/rrn)2.  so  that  ^(n,.Z.3)it  a;  104.  Therefore 

n,  2  _  104 

n  ~  (nL3)k 

Using  equation  (49) 

n^j(nlY.V.i*3xlO»a  =  2xlO'^ 

Define  the  [lacking  fraction  of  the  turbulence,  with  the  typical  volume  V3  = 
D,xPf/103rm3  of  beam-plasma  interaction. 

I  =  Jh-  - 

•V . .  10'Hi 

The  maxium  number  of  emitting  dipoles  which  could  fit  in  Vfrpt.  ,VmjJ.  is. 

v  _  V'lpl  _  10I5V3 
'  mor  4/3 ttZ.3  4(tM3 

I  he  actual  number  of  emitters  is.  using  equation  (19). 

•  - - i*e 


The  parking  fraction  is. 


/  a;  0.8- 


i 


□Xi 

Dai 

a 

I 

I 

| 

B 


meaning,  however.  Benford  (1987)  suggests  that  collisions  and  energy  exchang 
among  cavitons  can  produce  an  exp(  —  E2)  form  as  well. 


We  can  see  from  equation  (50)  that  the  packing  fraction  f  derived  from 
the  microwave  power  gives  good  agreement  with  the  experimentally  observed 
packing  fraction  of  strong  turbulence.  Since  the  microwave  power  caculation  is 
based  on  the  W( E)  distribution,  this  gives  support  to  our  long-range  correlation 
model. 

Note  that  we  have  assumed  a  correlation  existing  over  a  three  dimensional 
field,  even  though  E  may  be  a  one  or  two  dimensional  field.  We  picture  par¬ 
ticles  correlated  by  a  group  of  waves  which  can  be  strongly  anisotropic,  so  the 
correlation  integral  of  equation  (34)  yields  an  L3  form.  If  this  were  not  so.  the 
integral  would  be  proportional  to  A 2DL  for  a  one  dimensional  E.  and  A pL-  for 
a  two  dimensional  one.  Which  form  is  appropriate  for  a  given  experiment  must 
be  decided  by  details  of  the  case.  Clearly,  if  only  one  dimensional  correlations 
occur,  the  product  A 3DLnc  will  still  be  the  same  and  L  will  necessarily  be  much 
larger,  then  we  have  estimated  here. 
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Appendix  B 


Fluctuating  Electric  Fields  in  an  Intense 
Relativistic  Electron  Beam-Driven  Turbulent  Plasma 

Electrostatic  plasma  oscillation  in  plasmas  which  contain 
atoms  should  give  rise  to  satellite  lines  in  the  atomic 
spectrum.  The  satellite  lines  appear  near  the  forbidden 
transition. I ^ 

If  the  plasma  fields  were  time  independent  (or  with  rather 
low  frequency)  so  that  we  can  use  quasi-static  theory,  they  would 
produce  a  forbidden  line  through  the  static  Stark  effect  in  the 
well-known  manner,  i.e.,  the  two  satellites  coalesce  into  one 
forbidden  line.  Using  the  ratio  of  the  intensity  of  the 
forbidden  plus  satellite  lines,  to  the  allowed  line  intensity,  we 
could  derive  the  R.M.S.  field  as  the  combined  field  of 
oscillation  at  various  frequencies. 

The  level  system  in  Fig.  Bl,  according  to  [1],  has  usually 
different  levels  for  the  upper  states  of  both  lines.  Kawasaki  et 
al.,  ^2]  compare  forbidden  to  allowed  lines  which  have  tne  same 
initial  state. 

The  two  characteristics  shown  in  Fig.  B2  are: 

1)  We  need  not  make  any  assumption  (or  measurement!)  of  the 
population  density  of  the  forbidden  level,  since  the  lines  are 
emitted  from  the  same  level. 

2)  The  population  density  of  the  upper  level,  l,  can  be 
enhanced  by  a  resonance  absorption  from  the  level  A  .  For  making 
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the  resonance  absorption  effectively,  the  metastable  state  of  Hel 
is  suitable  for  the  level  i .  (Fig.  B3) 

The  application  of  tunable  dye  lasers  to  this  experiment 
offers  several  advantages : ^ 

1)  Spectroscopic  observations  integrate  along  the  line  of 
sight.  Satellite  intensities  are  always  measured  relative  to  the 
intensity  of  the  allowed  transition.  This  will  cause  serious 
errors,  especially  if  the  allowed  transition  is  emitted  from  the 
total  plasma  along  the  line  of  sight.  The  plasma  satellites, 
however,  stem  from  only  the  turbulent  region,  i.e.,  most  likely 
from  the  current-carrying  part  of  the  plasma.  The  laser- 
fluorescence  technique  enables  us  to  gain  spatial  resolution. 

(See  Fig.  B4) 

2)  The  laser  is  focused  into  the  plasma  and  resonantly 
tuned  to  the  allowed  transition.  This  improves  the  signal-to- 
noise  ratio,  as  well  as  the  time  resolution  of  the  measurement. 
Possible  impurity  lines  as  well  as  lines  from  molecules  can  be 
discriminated  against  by  pumping  the  allowed  transition  with  a 
narrow-band  laser  beam. 


In  this  way  we  have  succeeded  making  measurements  of  the 

o  1/2 


kV 


10  —  with  a  resolution  of 
cm 


electric  field  and  find  <E  > 
about  (mm) ’ . 

We  intend  to  go  on  measuring  the  E-field  in  different 
regions  in  the  plasma  —  on-axis  and  off-axis,  close  to  the 
anode,  and  downstream.  Also,  we  shall  measure  the  fields  in 
different  times  relative  to  the  electron-beam  pulse.  This  will 
yield  E(r,z,t).  We  expect  to  see  growth  and  dispersion  of 
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E(r,z,t),  which  will  yield  understanding  of  nonlnear  turbulence 
characteristics. 

A  sensitive  probing  of  this  low-E  region  can  make  contact 
with  conventional  "strong"  turbulence  theory  (0.1  <  W  <  1). 
Similarly,  beam-induced  fluorescence  will  yield  about  1000  times 
more  photons  at  high-E  Stark  shifts,  and  extend  our  measurement 
of  N  ( E )  to  300  kV/cm ,  with  the  same  statistical  reliability  as 
our  present  highest  measurement,  120  kV/cm.  This  will  provide 
check  for  any  break  in  the  exp(-E2)  distribution  (Appendices  A 
and  B) ,  and  thus  for  theoretical  models  of  high-E  turbulence. 

Rapid  switching  of  beam  current  and  voltage  will  allow  study 
of  high-E  lifetimes  in  the  plasma.  Our  preliminary  work 
(described  under  "Turbulence  Lifetime  Measurements"  in  our  former 
proposal)  shows  that  turbulence  persists  long  after  usual  theory 
predicts  its  dissipation.  We  plan  extensive  experiments 
switching  beam  current  on  to  explore  how  rapidly  high-E  regions 
appear,  and  similar  switch-off  experiments  to  measure  lifetimes. 

We  will  develop  theoretical  models  to  explain  the  wide  range 
of  experiments  we  have  already  done  or  will  soon  perform.  This 
will  include  a  careful  consideration  of  microwave  emission  data, 
and  comparison  with  turbulent  electromagnetic  emission  models. 
Coupling  of  the  microwave  measures  with  direct  N ( E )  data  will 
provide  a  valuable  constraint  on  theory. 


were 


3*  r 


a 


*»!l 


»>l 

'tj 

>?3 

I 


References : 


M.  Baranger  and  B.  Mozer,  "Light  as  a  Plasma  Probe."  Phy, 
Rev.  123 ,  25  (1961) 


K.  Kawasaki,  T.  Usui  and  T.  Oda ,  "Forbidden  Transitions  in 
Helium  and  Lithium  Due  to  Fluctuating  Electric  Fields  for 
Plasma  Diagnostics."  Journal  of  the  Phys.  Soc .  of  Japan,  51 
3666  (1982)  . 


H.J.  Kunze,  "Investigation  of  Plasma-Satellites  by  Laser- 
Fluorescence  Spectroscopy"  -  PJ.  517  in  "Spectral  Line 
Shapes"  ed .  B.  Wende  (1980). 


Au-owtQ 


Fo*ftiPO£A) 


,  M _ — i _ . 

Figure  B2 j  The  forbidden  line  is  caused  by  the 

level  mixing  of  i  with  j  due  to  electric 
fields,  while  the  allowed  line  is  due 
to  the  transition  from  i  to  1.  An  in¬ 
crease  in  population  density  of  i  enhances 
both  of  the  lines. 


orncAL 

tV/Af'UC,. 


f  F^ftBiOOEKJTIiAAJS  iTlOtJ 

4US2.X 


